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Abstract: We use the GEANT4 Monte Carlo framework to calculate the 7-ray albedo of the Moon due 
to interactions of cosmic ray (CR) nuclei with moon rock. Our calculation of the albedo spectrum agrees 
with the EGRET data. We show that the spectrum of 7-rays from the Moon is very steep with an effective 
cutoff around 4 GeV (600 MeV for the inner part of the Moon disc). Apart from other astrophysical 
sources, the albedo spectrum of the Moon is well understood, including its absolute normalisation; this 
makes it a useful "standard candle" for 7-ray telescopes, such as the forthcoming Gamma Ray Large Area 
Space Telescope (GLAST). The albedo flux depends on the incident CR spectrum which changes over the 
solar cycle. Therefore, it is possible to monitor the CR spectrum using the albedo 7-ray flux. Simultane- 
ous measurements of CR proton and helium spectra by the Payload for Antimatter Matter Exploration and 
Light-nuclei Astrophysics (PAMELA), and observations of the albedo 7-rays by the GLAST Large Area 
Telescope (LAT), can be used to test the model predictions and will enable the GLAST LAT to monitor 
the CR spectrum near the Earth beyond the lifetime of PAMELA. 



Introduction 

Interactions of Galactic CR nuclei with the atmo- 
spheres of the Earth and the Sun produce albedo 
7-rays due to the decay of secondary neutral pi- 
ons and kaons (e.g., [10, 7]). Similarly, the Moon 
emits 7-rays due to CR interactions with its sur- 
face [6, 11]; low energy 7-ray spectroscopy data 
acquired by the Lunar Prospector were used to 
map the elemental composition of the Moon sur- 
face [4, 8]. However, contrary to the CR interac- 
tion with the gaseous atmospheres of the Earth and 
the Sun, the Moon surface is solid, consisting of 
rock, making its albedo spectrum unique. 

Due to the kinematics of the collision, the sec- 
ondary particle cascade from CR particles hitting 
the Moon surface at small zenith angles develops 
deep into the rock making it difficult for 7-rays to 
get out. A small fraction of all produced pions, 
splash albedo pions, are mostly low energy ones 
thus producing the soft spectrum 7-rays. High- 
energy 7-rays can be produced by CR particles hit- 
ting the Moon surface in a close-to tangential di- 
rection. However, since it is a solid target, only the 



very thin limb contributes to the high energy emis- 
sion. 

The 7-ray albedo of the Moon has been calculated 
by Morris [6] using a Monte Carlo code for cas- 
cade development in the Earth's atmosphere that 
was modified for the Moon conditions. However, 
the CR spectra used as input in [6] differ con- 
siderably from recent measurements by AMS and 
BESS. Besides, due to the lack of accelerator data 
and models a number of approximations and ad- 
hoc assumptions were required to calculate the 
hadronic cascade development in the solid target. 

The Moon has been detected by the EGRET as a 
point source with integral flux F(>100 MeV) — 
(4.7 ± 0.7) X 10-"^ cm-2 g-i [11], -24% below 
the predictions [6]. The observed spectrum is steep 
and yields only the upper limit —5.7 x 10^^^ cm~-^ 
s^^ above 1 GeV. 

We report preliminary results for calculations 
of the 7-ray albedo from the Moon using the 
GEANT4 [1] framework code and discuss the con- 
sequences of its measurement by the upcoming 
GLAST mission. 
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Figure 1: Beam/target/detector setup for simulat- 
ing CR interactions in moon rock. The primary 
beam enters the moon rock target with incident po- 
lar angle 9p. Secondary 7-rays are emitted with 
polar angle 9. The detection volume surrounds the 
target. 

Monte Carlo simulations 

In the present work, we use version 8.2.0 of 
the GEANT4 toolkit. Figure 1 illustrates our 
beam/target/detector setup for simulating CR inter- 
actions in the Moon. The primary CR beam (pro- 
tons, helium nuclei) is injected at different incident 
angles into a moon rock target. We take the compo- 
sition of the moon rock to be 45% Si02, 22% FeO, 
11% CaO, 10% AI2O3, 9% MgO, and 3% TiOs by 
weight, consistent with mare basalt meteorites and 
Apollo 12 and 15 basalts [4, 2, 8]. A thin hemi- 
spherical detector volume surrounding the target is 
used to record the secondary 7-ray angular and en- 
ergy distributions in the simulation. 

The 7-ray yield dY^ {Ep, cos 9p)/dE^d cos 9 is 
calculated using the GEANT4 beam/target setup 
with a Monte Carlo method; here Ep is the kinetic 
energy per nucleon of the incident particle, 9p is 
the incident polar angle, E^ is the energy of 7-rays, 
and 9 is the polar photon emission angle. 

Figure 2 shows the secondary 7-ray yields inte- 
grated over all emission angles outward from the 
Moon surface for protons with Ep — 500 MeV and 
5000 MeV at incident angles cos 9p — 0.1 and 1, 
respectively. The shapes of the yield curves for dif- 
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Figure 2: 7-ray yield per proton interaction inte- 
grated over all emission angles from the Moon sur- 
face. Line-styles: red-dashed, cos9p — 1; blue- 
solid, cos^p = 0.1. Line-sets: lower, Ep — 500 
MeV; upper, Ep = 10000 MeV. 

ferent incident angles are very similar to each other 
for the case of low energy protons where the sec- 
ondary particles (pions, kaons, neutrons, etc.) are 
produced nearly at rest. In this case the 7-ray emis- 
sion is produced in a number of processes: pion 
and kaon decay, secondary electron and positron 
bremsstrahlung, and so forth. A considerable flux 
of 7-rays is produced in nuclear reactions such as 
neutron capture and nonelastic scattering [4, 8]; 
the features below ^10 MeV are due to nuclear 
de-excitation lines, where the most prominent con- 
tribution comes from oxygen nuclei. In the high 
energy case, the secondary distribution for pro- 
tons incident near zenith has a cutoff above ^^500 
MeV. Further away from zenith higher yields of 
secondary 7-rays are produced while the spectrum 
of 7-rays becomes progressively harder. This is 
a result of the cascade developing mostly in the 
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forward direction: for near zenith primaries, most 
high energy secondary 7-rays will be absorbed in 
the target, while a small fraction of produced pi- 
ons and kaons, splash albedo particles, mostly low 
energy ones, produce the soft spectrum 7-rays; fur- 
ther from zenith, the high energy secondary 7-rays 
will shower out of the Moon surface. 

Calculations 

To calculate the Moon albedo at an arbitrary solar 
modulation level, we use the local interstellar spec- 
tra (LIS) of CR protons and helium as fitted to the 
numerical results of GALPROP propagation model 
(reacceleration and plain diffusion models. Table 1 
in [9]); the CR particle flux at an arbitrary phase of 
solar activity at 1 AU can then be estimated using 
the force-field approximation [3]. 

To fit the LIS CR spectra we choose a function of 
the form: 

^ = JoVa.(i?, + 6.)--, (1) 

2—1 

where the flux units are mT^ s^^ sr^^ 
(GeV/nucleon)^^ and the parameter values 
are given in Table 1 . The latter are not unique and 
other sets could produce similar quality fits, but 
this does not affect the final results. 

Figure 3 shows the calculated total 7-ray albedo 
spectrum for CR protons and helium compared to 
the EGRET data for periods of lower solar (up- 
per solid: $ = 500 MV) and higher solar activity 
(lower solid: $ = 1500 MV). Taking into account 
that the exact CR spectra during the EGRET obser- 
vations are unknown, the agreement with the data 
is remarkable. The broken lines show the spectra 
from the limb (outer 5') and the central part of the 
disc (20' across) for the case of higher solar ac- 
tivity. As expected, the spectra from the limb and 
the central part are similar at lower energies (<10 
MeV); at high energies the central part exhibits a 
softer spectrum so that virtually all photons above 
^600 MeV are emitted by the limb. 

The clear narrow line at 67.5 MeV in Figure 3 is 
due to neutral pions decaying near rest. The line 
perhaps always exists due to the splash albedo par- 
ticles, but in case of a thin and/or gaseous target 
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Figure 3: Calculated 7-ray albedo spectrum of the 
Moon. Line-styles: black-solid, total; blue-dotted, 
limb - outer 5'; red-dashed, centre - inner 20'. Up- 
per solid line: $0 = 500 MV; lower solid line: 
$0 = 1500 MV. Data points from the EGRET 
[11] with upper and lower symbols corresponding 
to periods of lower and higher solar activity, re- 
spectively. The differential 1 year sensitivity of the 
LAT is shown as the hatched region. 

(usual in astrophysics) it is indistinguishable from 
the background 7-rays. For the case of a thick solid 
rock target, and thus fast energy dissipation of the 
particle cascade, and because of the kinematics of 
the interaction, the continuum 7-ray background is 
much smaller revealing the narrow line from pion 
decay at rest. 

Discussion and Conclusions 

The GLAST LAT is scheduled for launch by 
NASA in winter of 2008. It will have superior an- 
gular resolution and effective area, and its field of 
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Table 1: Fits to local interstellar CR spectra 



view (FOV) will far exceed that of its predecessor 
EGRET [5]. The LAT will scan the sky continu- 
ously providing complete sky coverage every two 
orbits (--3 hr). About 20% of the time the Moon 
will be in the FOV at different viewing angles. The 
point spread function (PSF) of the instrument is 
~0.8°, at 1 GeV, but reduces dramatically at higher 
energies: ~0.5° at 2 GeV and ~0.2° at 10 GeV. 

The Moon with its steep albedo spectrum presents 
almost a black spot on the 7-ray sky above ^4 
GeV. The central part of the Moon has an even 
steeper spectrum with a cutoff at ~600 MeV. The 
albedo spectrum of the Moon is well understood 
while the Moon itself is a "moving target" passing 
through high Galactic latitudes and the Galactic 
centre region. This makes it a useful "standard can- 
dle" for the GLAST LAT at energies below 1 GeV. 
A simultaneous presence of the RAMELA on-orbit 
capable of measuring protons and light nuclei with 
high precision provides a necessary input for ac- 
curate prediction of the albedo flux and a possible 
independent calibration of the GLAST LAT. An 
additional bonus of such a cahbration is the pos- 
sibiUty to use GLAST observations of the Moon to 
monitor the CR spectra near the Earth beyond the 
projected 3 yr hfetime of the RAMELA. 

The line feature around 67.5 MeV from 7r°-decay 
produced by CR particles in the solid rock target 
is interesting. The lower energy Umit of the LAT 
instrument is below 20 MeV while the energy reso- 
lution is ~15% at 100 MeV, and improves at higher 
energies. With a suitable event selection it may be 
possible to observe the Une. There is no other as- 
trophysical object predicted to produce such a nar- 
row line and there is no other Une expected except, 
perhaps, from dark matter annihilation. A pos- 
sibility for energy calibration at higher energy is 
provided by the steep albedo spectrum above 100 
MeV: a small error in the energy determination will 
result in a large error in the intensity. 
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